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ABSTRACT 


A  method  is  proposed  for  whitening  the  spectra  of 
teleseismic  P  waves  to  aid  in  the  determination  of  event 
depths.  The  method  uses  a  network  of  common  recording 
stations,  all  of  which  must  have  recorded  signals  from 
a  suite  of  closely-spaced  events,  e.g.,  a  set  of  under¬ 
ground  explosions.  Individual  spectra  are  whitened  by 
removing  a  path-recording  site  response  determined  by 
averaging  spectra  from  the  suite  of  events  as  recorded 
at  a  single  station.  The  whitened  spectra  from  a  single 
event  are  then  summed  over  the  recording  network  to 
yield  a  whitened  source  spectra.  The  whitening  process 
enhances  spectral  components  in  the  band  2  to  5  Hz. 
Application  of  the  proposed  whitening  method  to  explo¬ 
sion  spectra  for  theoretical  and  experimental  data  sets 
indi;.ates  that  the  method  may  be  useful  in  analyzing 
spectral  nulls  provisionally  interpreted  as  due  to  pP-P 
interference.  Depth-phase  analysis  for  seven  NTS  explo¬ 
sions,  using  both  whitened  and  unwhitened  data,  shows 
that  spectral  nulls  are  enhanced  in  the  whitened  spectra. 
Network  summed  spectra  (unwhitened)  for  seventeen  earth¬ 
quakes  are  also  included;  these  demonstrate  the  fact  that 
nulls  can  be  caused,  for  example,  by  multipath  effects  and 
multiplet  sources  as  well  as  by  depth-phase  interference. 
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INTRODUCTION 


of  this  work  is  to  demonstrate  a 

rilo  a  teleseismic  spectra.  As  applied  to 

explosion  data,  the  method  can  enhance  spectral  charac- 
teristics  associated  with  interference  of  the  type  pP-P 
and  thus  aid  in  source  depth  analysis.  Previous  depth-  ’ 
ermination  analyses  (Cohen,  1970)  have  employed 
in  ividual  station  spectra,  which  in  an  attempt  to 
enhance  spectral  components  in  the  band  2  to  5  Hz  were 
corrected  for  instrument  response.  Further,  to  reduce 
reverberation  noise  at  the  individual  recording  sites 

DeslTtMs'  “  "^■•'-'’'-oraged  spectrum.’ 

Despite  this  processing,  the  network-summed  spectra 

were  still  band-limited  (primarily  by  the  earth’s 

TZZTr  generally  exhibited 

e  single  spectral  null  which  could  be  interpreted  as 

deriving  from  depth-phase  interference. 

i-t  spectral  characteristics  in  the  band 

z,  eleseismic  spectra  requite  whitening.  Here  we 
whiten  the  spectra  by  dividing  out  a  path-recording 
ite  response  which  is  determined  by  averaging  spectra 
from  a  suite  of  events  as  determined  from  data  recorded 
a  single  station.  The  whitened  spectra  from  a  single 
event  are  then  summed  over  the  recording  network  to 
y  e  the  average  spectra  u  ed  in  depth-phase  analysis. 

snect!^'^"'"  technique  can  enhance, 

spectral  characteristics,  we  must  add  that  scalloped 

spectra  can  be  caused  by  multi-path  effects  and  multi- 

P  e  sources  as  well  as  by  depth-phase  interference. 
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Examples  of  scalloped  earthquake  spectra  are  shown  in  the 
Appendix.  Thus,  caution,  must  be  exercised  in  using 
spectral-null  depth  determinations  as  a  diagnostic. 
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METHOD 


Let 

E.(f)  =  source  spectrum 

S^Cf)  =  explosion  site  crustal  response 
P,  . (f)  =  mantle  path  response 
R.  (f)  =  receiving  site  crustal  response 
I . (f)  =  instrument  response 

Then,  for  the  i'th  underground  shot  as  recorded  at  the 
j'th  station,  the  observed  short-period  spectrum  Cf) 

is  given  by: 

-  [1  +  +  20005(^X3^3)] » 

where  the  quantity  in  square  brackets  is  the  pP-P  inter 
ference  function,  o  is  the  free-surface  reflection  coef 
ficient,  and  is  the  depth-phase  delay  time  at  the 
j'th  station  for  the  i'th  shot.  Note  that: 

t;..  ~  and  can  be  taken  simply  as 

however, 


since  in  this  work  we  will  use  a  suite  of  shots  having 
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various  source  depths.  In  previous  work  Cohen  (1970) 
used  a  network-average  spectrum  to  enhance  pP-P  spec¬ 
tral  periodicities.  Further,  the  individual  station 
spectra  were  corrected  for  instrument  response.  The 
network-average  spectrum  used  for  depth-determination 
analysis,  therefore,  is  given  by: 

?i(f)  =  ^ 

2 

•  [1  +  a  +2  cos(o)T^)] 


(2) 


where  the  sum  is  performed  over  the  spectra  at  N  record¬ 
ing  sites.  Implicit  in  the  above  relation  is  the  assump¬ 
tion  that  'holes'  in  the  spectra  which  are  caused  by 
the  geology  at  the  individual  recording  sites  and  along 
the  source-receiver  paths  will  in  large  part  average  out. 
If  this  assumption  holds,  and  both  E^(f)  and  S^(f)  change 
slowly  over  the  band  0,3-5  Hz,  the  character  of  <j)j|^(f) 
will  be  controlled  by  the  cosinusoidal  depth-phase 
function.  That  is: 


^i(f) 


»  E^(f)  ‘S^Cf)  •[!  +  + 


2acos  (u)T^)  ]  ‘TO, 


(3) 


where 


Though  application  of  this  technique  has  met  with 
some  success  (Lambert  et  al.,  1969;  Cohen,  1970),  the 
method  is  not  withovt  its  shortcomings.  While  we  may 
be  reasonably  assured  that  E^(f)  varies  slowly  in  the 
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frequency  band  of  interest  (cf.  Wyss  et  al.,  1971), 

Wirth  (1971)  has  shown  that  for  interval  reflection 
coefficients  somewhat  more  extreme  than  are  usually 
found,  teleseismic  spectra  can  be  considerably  compli¬ 
cated  due  to  the  effects  of  S^(f).  Further,  the  effect 
of  (f)  for  every  shot-station  pair  is  approximately 
the  same,  generally  to  attenuate  the  higher  frequencies; 
the  result  is  that  short-period  spectra  are  severely 
band-limited.  Finally,  from  a  practical  point  of  view, 
a  large  number  of  recording  stations  --  say  10  --  is 
not  always  available  for  network  averaging,  and  so  the 
^^(f)  may  still  reflect  contributions  from  the  Rj(f), 

In  sum,  the  network-averaged  spectrum  7j^(f)  tends  to  be 
complicated,  showing  predominance  of  energy  in  the 
0.5  -  2.0  Hz  band,  and  often  exhibits  only  a  single 
spectral  null  which  can  be  attributed  to  interference 
between  P  and  pP.  Two  examples  supporting  the  above 
observations  are  given  in  Figures  la  and  lb.  In  the 
case  of  GASBUGGY,  the  P^^  spectra  for  12  stations  were 
averaged;  for  LONG  SHOT,  31  station  spectra  were  used. 

In  both  cases  the  network -averaged  spectra  are  band- 
limited  and  exhibit  a  single  broad  spectral  null  which 
has  been  interpreted  as  the  depth-phase  null. 

Suppose  now  that  we  have  a  network  consisting  of 
several  seismic  stations,  each  of  which  records  the 
signals  from  a  group  of  explosions  in  the  same  general 
epicentral  region.  Further,  let  the  explosions  be  at 
various  depths.  Summing  over  spectra  at  a  single  station, 
we  have: 

♦  .(f)  -  ^  {J  E.(f).S.(f)-P.  .(f)-[l  ♦  ♦  2a  cos(u)t.  .)  ]  }*R.  (f) 

(4) 
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when  the  sum  is  performed  over  M  explosions. 
But 


♦  .(f)  =  ^  E(f)-P.(f)^R^.(f) 

tlS.  (f)[l  +  ♦  2ocos(u)T  .  .)  ] } 

i  1  ij 

since  in  the  frequency  band  of  interest, 

E^(f)=Ejjj(f)  -  E(f),  and  the  path  from  the  explosions  to 
a  given  receiver  are  virtually  identical  (i.e,, 

Pij  (f)sPmj (f)  “  Pj(f))«  With  the  exception  of  spectral 
characteristics  due  to  gross  crustal  structures  at  the 
source, we  would  expect  the  factor 


(5) 


^  JS^(f)»[l  ♦  ♦  2ocos(wTjj)l  E  SU  (6) 

to  be  a  smoothly  varying  function. 

Thus 


♦  .(f)  -  E(f)«Pj(f)*Rj(f).SlT  (7) 

Let  us  now  use  the  ♦j (f)  to  correct  the  raw 
observed  spectra.  Again,  instrument  response  has  been 
removed,  although  it  would  have  canceled  anyway  in  the 
case  of  identically  equipped  stations.  Thus: 

♦  E(f).S.(f).P,(f).R,(f).(l  .  *  2ocos(ut,J) 

e_(f)-.JJ - .  - i - 2 - 2 - _i2 _ 

E(f)-'p.(f)-Rj'(f)-SIT  (3) 

■  ~  S.  (f)  (1  ♦  ♦  2ocos(wt..)] 

STT  ^ 
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Averaging  the  e^jCf)  over  the  network  of  recording  sta 
tions  yields 


7i(f) 


— ♦  nacosCwT.  .)] 
N-STT  ^  j 


S^(f) 

“  -^TT” 


[1  +  a' 


2acos(a)T^)  ] , 


where  (S.(f)/SI5)  should  now  contain  only  the  more  minor 
contribuiions  from  S..  that  is,  the  contributions  due  to 
small-scale  layered  structures  particular  to  a  given 
explosion  site# 

Comparing  equation  (9)  with  equation  (3),  we  see 
that  the  proposed  analysis  method,  through  the  removal 
of  source  and  path  responses,  may  be  expected  to  yield 
a  whitened  network-average  spectrum,  enriched,  for 
example,  in  the  scalloped  components  produced  by  p 
interference. 


(9) 
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TEST  OF  THE  PROPOSED  METHOD 


We  first  test  the  proposed  analysis  method  on 
simulated  teleseismic  data.  Specifically,  we  will  gen¬ 
erate  short-period  seismograms  at  seven  hypothetical 
receiving  stations  for  seven  hypothetical  underground 
explosions.  The  receiving  stations  are  assumed  to  have 
a  world-wide  distribution,  and  so,  the  crustal  sections 
from  station  to  station  can  exhibit  large-scale  differ¬ 
ences.  On  the  other  hand,  the  explosions  are  assumed 
detonated  in  a  region  of  small  dimensional  extent  and 
as  such,  only  very  minor  changes  in  velocities  and 
layer  thicki.ess  are  permitted  from  site  to  site. 

The  model  used  to  simulate  a  nuclear  explosion 
has  a  power  spectrum  given  by  Haskell  (1967).  Using 
this  model,  a  minimum  phase  signal  for  a  100  kT  explo¬ 
sion  in  granite  was  computed.  This  signal  was  then 
filtered  to  incorporate  the  short-period  Benioff  re¬ 
sponse  (pendulum  and  glavonmeter  damping  taken  as  1.0; 
pendulum  free  period:  1,25  seej  galvonoiaeter  free  period 
0,2  sec).  The  resulting  wavelet,  which  was  convolved 
with  the  spike  seismograms  corresponding  to  the  various 
crustal  models  employed,  is  shown  in  Figure  2.  A  Monte 
Carlo  technique  was  use^  to  generate  layer  models  for 
both  the  source  and  receiver  sites.  Basic  to  each 
crustal  model  was  the  realistic  gross  velocity-depth 
distribution  shown  in  Figure  3.  Random  variations  to 
this  model  were  introduced  to  produce  the  crustal 
structures  at  the  source  and  receiver  sites  (Figures  4, 
and  5  through  8,  respectively").  Because  the  explosions 
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«cre  assumed  detonated  in  the  same  general  area,  only 
,,inor  variations  in  the  velocities  and  layer  thickness 
of  a  single  gross  crustal  model  were  “ 

source  site  to  another.  As  such,  the  crustal  models  for 

the  shot  sites  are  virtually  identical  to 

Figure  4.  On  the  other  hand,  the  crustal  models  for  the 

seven  receiver  sites  display  marked  differences, 

is  expected  for  sites  located  in  widely  separated 

geological  provinces. 

Given  the  layered  models,  synthetic  seismograms 
were  computed  which  took  into  account  all 
occurring  at  interfaces  to  approximately  200  km  depth 
:r::th  the  source  and  receiver.  An 
for  the  combined  source  and  receiver  stac  gi 
Wirth  (1971)  was  used  to  generate  the  synthetic  seismo¬ 
grams.  The  program  generates  solutions  for  a  spe c  d 

number  of  receivers  in  fixed  positions  for  a 
number  of  explosion  sources  of  given  dep  hs  Depths 
the  hypothetical  events  used  here  are  listed  in  Table  1, 
to  ether  with  the  frequency  of  the  first  spectra  nu  _ 
for  a  nominal  overburden  velocity  of  5  ^  J 

thetic  seismograms  so  generated  are  shown  in  Figure 

through  15. 

Application  of  the  depth-analysis  technique 
described  by  Cohen  (1970)  and  that  proposed  here  yiel 
the  results  shown  in  Figures  16  through  19.  Spectra 
have  been  smoothed  with  a  15-point  running  average 
enhance  the  display  of  broad  nulls  which  might  b 

associated  with  pP.  Frequencies  where 

ore  summarized  in  Table  II,  together  with  implied 
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depths  based  on  a  nominal  over-burden  velocity  of 
3  km/sec.  Actually,  average  overburden  velocities  for 
the  hypothetical  shots  varied  from  3,1  to  3.5  km/sec, 
depending  on  the  shot  depth.  As  such,  depth  estimates 
in  Tables  I  and  II  are  somewhat  shallow. 

As  seen  from  Table  II,  the  depth-phase  analysis 
method  employing  unwhitened  data  yielded  results  in 
only  three  cases;  in  each  of  these  cases,  the  depth  of 
the  event  is  greater  than  1.0  km.  The  proposed  depth 
analysis  method,  however,  yielded  satisfactory  results 
in  every  case,  with  shot  depths  ranging  from  0,7  to 
1.4  km.  It  must  be  admitted,  however,  that  a  plot  of 
the  unwhitened  spectra  on  a  log  scale  would  al'^o  show 
at  least  a  significant  first  null  in  every  cas2.  k'e 
see  that  defining  the  word  "significant"  is  difficult. 
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RESULTS  -  experimental 

Application  of  the  proposed  analysis  method  to 
real  data  presents  some  problems.  In  the  previous 
relrtlical  example,  every  station  ..recorded"  signals 
from  every  shot;  we  find  for  real  data  that  the  more 
stations  we  consider,  the  fewer  common  shots  we 
available  (and  vice  versa) .  This  of  “ 

consequence  of  not  having  some  given  set  of  LRSM  an 
observatory  stations  continuously 

,  =  it  is  still  possible  to  test  the  method  using 

h  ;a  C  ed  at  seven  stations  for  three  to  seven 
Z  eventrstation  locations  relative  to  NTS  are  given 
in  Figure  20,  while  pertinent  shot  data  is  liste  in 
rlble  IV  lists  the  Shot  data  used  for  each 

Station. 

inspection  of  Table  III  shows  that  four  of  the 

S  rs  to  be  used  have  depths  in  the  range  .55  to 
seven  shots  to  oossibility  of  weighting 

67  km  Thus,  we  considered  the  possiDix  y 

Z  s;;ctra  or  these  four  shots  such  that  the  spectra 
the  spectra 

sum  over  the  sui^e  ot  even 

net  be  biased  by  the  characteristics  of  ^^ots  in  this 
depth  range.  That  is.  we  were  somewhat  apprehens 
about  the  effectiveness  with  which  the  cosinusoi 
depth  phase  functions  would  average  ‘ 

-  -  :r:;-  - l;::: 

:s;;:.“r::rh:'we'decided .  weight  the  shot  spec, 
tra  equally. 

Average  station  spectra  obtained  by  summing  over 
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the  available  shot  spectra  at  a  given  station  are  shown 
in  Figure  21,  Except  for  the  average  spectra  at  CPO  and 
KC-MO,  the  average  station  spectra  show  a  single,  broad 
spectral  rise  with  peak  power  observed  in  the  band 
,7  -  ,9  Hz,  The  spectra  at  three  stations  --  HN-ME, 
RK-ON,  and  SV3QB  --  exhibit  slight  nulling  between 
,9  -  1,0  Hz,  but  these  nulls  are  considered  second-order 
characteristics.  Presumably  then,  the  average  station 
spectra  ♦j(f)  at  HN-ME,  PG-BC,  RK-ON,  SV3QB,  and  WMO 
primarily  exhibit  the  product  of  path  and  recording  site 
crustal  responses  for  each  station  (equation  (7)), 

The  average  station  spectra  at  CPO  and  KC-MO  are 
somewhat  more  complicated  than  the  spectra  at  other 
stations.  Whether  we  are  seeing  real  perturbations 
caused  by  path  and  site  responses,  or  merely  the 
results  of  averaging  over  an  insufficient  number  of 
shots  (CPO: 3;  KC-MO: 5)  is  a  matter  of  conjecture. 

Because  of  the  shot-depth  distributions  for  events 
recorded  at  CPO  and  KC-MO  (Table  III),  we  would 
expect  the  average  station  spectra  at  these  stations 
to  exhibit  some  nulling  effects  due  to  shots  with 
depths  between  0,55  and  0,67  km.  However,  we  will 
proceed  with  the  analysis  using  all  of  the  spectra 
shown  in  Figure  21,  again  emphasizing  our  blindness  to 
shot  depth  distribution. 

Spectral  null  analyses  performed  on  the  raw  spectra 
and  on  those  spectra  corrected  for  path  and  site 
response,  are  shown  in  Figures  22  through  35,  The 
results  are  presented  in  pairs,  grouped  by  shot.  The 
bottom  trace  in  each  case  is  the  average  network 
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spectrum  computed  over  the  suite  of  individual  station 
spectra. 

BOXCAR  (Figures  22  and  23) 

A  conventional  null- frequency  analysis  (Cohen,  1970) 
performed  using  spectra  at  five  LRSM  stations  is  shown 
in  Figure  22.  All  station  spectra  show  a  sharp  drop  in 
power  at  .9  Hz,  with  nulls  clearly  visible  near  this 
frequency  in  che  spectra  at  HN-ME  and  SV3QB.  A  sum 
over  the  station  network  exhibits  a  null  at  .93  Hz. 
Provisionally  interpreting  this  null  to  derive  from 
pP-P  interference,  and  taking  .93  Hz  as  the  null- 
frequency  increment,  the  corresponding  depth-phase 
delay  time  is  1.07  seconds. 

Analysis  of  the  station  spectra  corrected  for 
path  and  site  response  (Figure  23)  lends  credence  to 
the  above  result.  All  station  spectra  now  exhibit  a 
null  near  .9  -  1.0  Hz.  Further,  several  spectra  show 
a  second  drop  in  power  near  1.8  Hz,  with  a  null  at 
this  frequency  observed  in  the  SV3QB  spectrum.  The 
network- summed  spectrum  exhibits  a  null  at  0.95  and  a 
weak  null  at  1.80  Hz,  suggesting  a  null-frequency 
increment  of  0.9  Hz.  This  yields  an  estimated  depth- 
phase  delay  time  for  the  event  of  1.1  seconds. 

CHARTREUSE  (Figures  24  and  25) 

Conventional  depth-phase  analysis  using  the  spectra 
at  six  stations  is  shown  in  Figure  24.  All  spectra 
show  a  significant  drop  in  power  near  1.3  -  1.4  Hz, 
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with  nulls  observed  near  these  frequencies  at  KC-MO, 
PG-BC,  and  SV3QB,  The  network  sum  shows  a  slight 
depression  at  about  1.45  Hz,  Taking  1.45  Hz  as  the  null- 
frequency  increment  for  interference  of  the  type  pP-P 
yields  a  depth-phase  delay  time  of  .6t  second. 

Results  obtained  using  corrected  station  spectra 
(Figure  25)  are  somewhat  ambiguous.  While  broad  nulls 
are  observed  in  the  average  spectrum  at  1,4  and  2.6  Hz, 
one  might  argue  that  these  nulls  are  the  second  and 
fourth  in  a  series  have  a  null-frequency  increment  of 
,6  to  ,7  Hz,  However,  nulls  at  0,6  and  1.9  Hz  are 
poorly  developed.  These  data  then,  suggest  a  null- 
frequency  increment  of  1.3  Hz,  with  a  corresponding 
depth-phase  delay  time  of  ,77  second. 

DURYEA  (Figures  26  and  27) 

Conventional  spectra-null  analysis  yields  the 
results  shown  in  Figure  26,  A  broad  null  near  1,7  Hz 
is  observed  in  the  network -summed  spectrum,  though 
this  feature  is  not  well  developed.  As  such,  the  data 
suggest  a  null-frequency  interval  of  1,7  Hz,  with  a 
corresponding  delay  time  of  ,59  second. 

Null-frequency  analysis  for  DURYEA  using  spectra 
corrected  for  path  and  site  response  (Figure  27)  yields 
a  similar  result.  Again,  a  null  in  the  average  spectrum 
is  observed  near  1.7  Hz,  with  another  poorly  developed, 
broad  null  around  3,25  Hz,  Nulls  at  1,7  Hz  are  observed 
in  the  spectra  for  CPO,  mN-ME  and  SV3QB  with  the  spec¬ 
trum  at  the  latter  station  also  exhibiting  a  null 
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at  3.4  Hz.  For  this  shot,  then,  nulls  provisionally 
interpreted  as  due  to  pP-P  interference  are  observed 
with  possible  null-frequency  increments  of  1.7  Hz, 
yielding  a  corresponding  estimate  for  the  depth-phase 
delay  time  of  0.59  second. 


GREELEY  (Figures  28  and  29) 

Inspection  of  the  spectra  obtained  through  conven¬ 
tional  null -frequency  analysis  (Figure  28)  shows  several 
stations  to  exhibit  nulls  near  1  Hz.  The  null  at  1  Hz 
observed  at  SV3QB  is  particularly  well  developed.  Upon 
summing  over  station  spectra,  however,  we  obtain  an 
average  spectrum  devoid  of  nulls.  As  such,  spectral- 
null  analysis  using  the  unprocessed  station  spectra  is 

not  recommended. 

Recourse  to  the  processed  station  spectra  (Figure 
29)  yields  a  most  satisfactory  result.  A  spectral  null 
in  the  vicinity  of  1  Hz  is  observed  at  every  station, 
with  the  spectra  at  HN-ME  and  PG-BC  also  exhibiting  weak 
nulls  near  2  Hz.  We  are  therefore  inclined  to  adopt  a 
null-frequency  interval  of  1  Hz,  with  a  corresponding 
depth-phase  delay  time  of  one  second. 

HALF- BEAK  (Figures  30  and  31) 

Conventional  depth-phase  analysis  for  HALF-BEAK 
(Figure  30)  yields  a  questionable  result.  Spectra  at 
KC-MO  and  RK-ON  exhibit  weak  nulls  around  1.2  Hz, 
with  a  very  weak  null  observed  in  the  network- 
averaged  spectrum  at  1.25  Hz.  Though  speculative,  we 
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adopt  1.25  Hz  as  the  null- frequency  interval,  yielding 
a  depth-phase  delay  time  of  0.80  second. 

Using  spectra  corrected  for  'Sth  and  site  response 
(Figure  31)  yields  a  somewhat  more  reliable  estimate  of 
the  null -frequency  interval.  Here,  a  null  in  the  network- 
averaged  spectrum  is  clearly  observed  at  1.15  Hz,  with 
a  weak  null  present  around  2.0  Hz.  Taking  1.15  Hz,  as 
the  average  null-frequency  interval,  we  compute  a  depth- 
phase  delay  time  of  0.87  second. 

KNICKERBOCKER  (Figures  32  and  33) 

Conventional  depth-phase  analysis  (Figure  32) 
suggests  the  presence  of  a  broad  null  centered  at  1.4  Hz, 
perhaps  best  developed  in  the  spectrum  for  HN-ME.  As 
such,  the  null -frequency  interval  will  be  taken  as 
1.4  Hz,  yielding  a  depth-phase  delay  time  of  0,7  second. 

Employing  the  proposed  spectral  whitening  technique 
yields  an  interesting  result  (Figure  33).  Nulls  in  the 
network- averaged  spectrum  are  observed  at  1,30  and 
2.15  Hz  with  nulls  at  these  frequencies  best  developed 
in  the  station  spectra  for  CPO,  HN-ME,  PG-BC  and  WMO. 
Further,  the  CPO  spectrum  also  exhibits  nulls  at  0.40 
and  2.85  Hz.  This  suggests  that  the  null  frequencies  are 
given  by  the  formula: 

=  (2n+l)fQ,  n  =  0,1,2,,.., 

where  f  «  0.42  Hz.  But  this  implies  that  the  null 
phenomenon  observed  here  is  caused  by  the  interference 
of  two  similar  signals  arriving  with  the  same  first-motion 
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polarity,  with  the  second  signal  delayed  1/ (2£q)  or  1.2 
seconds  relative  to  the  first  signal.  The  second  signal 
presumably  is  produced  by  spalling  o£  the  near- surface 
layers  (Eisler  and  Chilton.  1964),  and  will  therefore 
be  referred  to  as  Pgpy^LL* 

In  the  case  of  KNICKERBOCKER,  then,  spectral  null¬ 
ing  is  apparently  controlled  by  interference  of  the  type 
P  „-P.  with  interference  of  the  type  pP-P  less  well 
developed,  if  developed  at  all.  As  such,  the  depth  esti¬ 
mate  previously  obtained  using  the  unwhitened  KNICKERBOCKER 
data  is  suspect,  for  the  null  at  1.4  l\z  observed  in  the 
network-averaged  spectrum  (Figure  32)  is  probably  related 

PSPALL-P  interference. 

REX  (Figures  34  and  35) 

The  network- average  spectrum  obtained  using  con¬ 
ventional  analysis  techniques  (Figure  34)  yields  an 
ambiguous  result.  A  well-developed  null  is  observed  at 
1.4  Hz,  and  taking  this  frequency  as  the  null  frequency 
interval  yields  a  depth-phase  delay  time  of  0.7  second., 

On  the  other  hand,  a  weak  null  in  the  average  spectrum 
is  observed  at  2.3  Hz,  with  the  suggestion  of  a  null 
also  seen  at  0.45  Hz.  These  nulls,  together  with  the 
null  at  1.4  Hz,  suggest  a  null  frequency  interval  of 
0.9  Hz,  with  nulling  presumably  due  to  PgpALL'P 
ference.  If  the  latter  explanation  is  correct,  Pgp^LL 
arrives  1.1  seconds  after  the  P  arrival. 

Application  of  the  spectral  whitening  technique 
yields  a  negative  result.  Numerous  nulls  are  observed 
in  the  individual  station  spectra,  but  little 


consistency  in  the  null  frequencies  is  observed  from 
station  to  station.  As  such,  it  is  not  surprising  to 
find  that  the  network-averaged  spectrum  displays  little 
character.  The  spectra  for  several  stations  (PG-BC, 
RK-ON,  and  SV3QB)  do  exhibit  broad  peaks  centered  near 
3.3  Hz,  with  addendant  nulling  suggested  in  the  band 
1.7-2  Hz.  Whether  this  broad  nulling  is  due  to  pP-P  or 
Pspall"P  interference,  however,  is  purely  speculative. 
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DISCUSSION 


Depth-phase  analysis  resuls  are  shown  in  Figure  36. 
Results  for  KNICKERBOCKER  and  REX  are  not  included.  To 
facilitate  a  determination  of  avi^rage  overburden  velocity, 
the  delay  times  are  plotted  against  double  the  known  shot 
depths. 

Both  the  whitened  and  un’yhitened  spectra  yield  a 
consistent  set  of  delay-time  determinations.  In  general, 
nulls  in  the  spectra  corrected  for  site  and  path  response 
are  more  enhanced  than  those  observed  in  spectra  not 
corrected  for  these  effects.  The  data  suggest  an  average 
overburden  velocity  of  1.8  km/sec  to  depths  on  the  order 
of  0.6  km,  with  an  increase  to  3  km/sec  suggested  at 
greater  depths.  This  results  in  good  agreement  with  the 
upper  crustal  model  for  the  Yucca  Flat  area  of  NTS  given 
by  Glover  and  Alexander  (1970).  Their  model  shows  a 
surface  layer  0.2  km  thick  with  a  compressional  velocity 
of  1.3  km/sec  followed  by  a  layer  0.8  km  thick  with  a 
velocity  of  3.0  km/sec. 
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CONCLUSIONS 


1,  Teleseismic  spectra  can  be  whitened  by  dividing 
out  a  path-recording  site  response  determined  by  averag¬ 
ing  spectra  from  a  suite  of  events  as  determined  from 
the  data  at  a  single  station.  Source  spectral  character¬ 
istics  can  then  be  enhanced  by  averaging  the  individual 
whitened  spectra  over  the  recording  network.  Implicit 
in  this  type  of  analysis  is  the  requirement  that  a  net¬ 
work  of  common  stations  has  recorded  signals  from  a 
suite  of  events. 

2.  The  whitening  process  tends  to  enhance  spect.al 
components  in  the  band  2  -  5  Hz.  As  such,  the  method 

may  be  useful  in  the  analysis  of  spectral  characteristics 
(nulls)  provisionally  interpreted  as  due  to  pP-P  inter¬ 
ference. 

3.  Analysis  of  short-period  data  for  seven  NTS 
explosions  shows  that  the  sought-after  depth-phase 
spectral  nulls  can  be  enhanced  through  application  of 
the  proposed  spectral  whitening  technique.  The  whitened 
spectra,  however,  are  somewhat  more  complicated  than 
the  unwhitened  spectra,  and  accordingly,  are  more 
difficult  to  interpret. 

4.  Recognizing  that  both  pP-P  and  PsPALL’^  inter¬ 
ference  can  be  present  in  teleseismic  P  spectra  suggests 
that  more  objective  techniques  for  examining  spectral 
characteristics  (e.g.,  the  pseudo-autocorrelation  and 
the  cepstrurc}  be  employed  in  null-frequency  dnalyris. 
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TABLE  I 


Shots  Depths  and  Nominal  First-Null  Frequenci 


SHOT 

DEPTH 

!L 

1 

.708 

2.12 

2 

1.349 

1.11 

3 

.772 

1.94 

4 

.883 

1.70 

5 

1.169 

1.28 

6 

.695 

2.16 

7 

1.101 

1.36 

^Overburden  velocity  of  3,0  km/sec  assumed 


TABLE  II 


Observed  Null  Frequencies  and  Implied  Depths 


SHOT  and 
METHOD* 

iL. 

Ji. 

1-A 

mm 

-- 

-- 

1"B 

2.50 

4.70 

-- 

2-A 

1.25 

-- 

-- 

2-B 

1.20 

2.45 

3.85 

3-A 

-- 

-- 

-- 

3-B 

1.95 

3.95 

-- 

4-A 

-- 

-- 

-- 

4"  B 

1.85 

3.80 

m  m 

S-A 

1.45 

•  • 

-- 

5-B 

1.40 

3.00 

4.20 

6-A 

-- 

m  m 

m  m 

6-B 

2.40 

(5.00) 

-- 

7-A 

1.45 

m  m 

-- 

7-B 

1.40 

2. 95 

4.15 

*A  -  Cohen  (1970) 

B  -  Present  Paper 

IMPLIED 

.  DEPTH  (km) 

*  IT  -  3  km/sec 
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DEPTH  (km) 
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0.71 

2.35 
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RECCBO  OF  SEISMOGRAM  AVAILABILITY 
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Figure  1.  Spectral  null  analysis  a)  GASBUGGY  b)  LONG  SHOT, 
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Figure  2.  Minimum-phase  wavelet. 
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Figure  4.  Crustal  model  for  the  source  sites 
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Figure  5.  Crustal  models  for  receiving  sites  1  and  2. 
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Figure  6.  Crustal  models  for  receiving  sites  3  and  4 
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Figure  8.  Crustal  model  for  receiving  site  7. 
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Figure  15.  Theoretical  station  scisaograas.  Shot 


Figure  15.  Theoretical  station  seismogram.  Shot 
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Figure  16.  Network-summed  spectra.  Shots  1  and  2 
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Figure  17.  Network-summed  spectra,  Shots  3  and  4. 
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Figure  18.  Network- summed  spectra,  Shots  5  and  6. 
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Figure  19.  Network -summed  sptictra,  Shot  7, 


gure  21.  Sum  spectra  at  individual  recording  stations 
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Figure  22.  Unwhitened  station  and  network-summed  spectra 
BOXCAR. 
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Figure  23.  Whitened  station  and  network  sunmed  spectra 
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Figure  24.  Unwhitened  stati 
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Figure  28.  Unwhitened 
GREELEY. 
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Figure  30.  Unwhitened  station  and  network-sumed  spectra  - 
HALF  BEAK. 
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Figure  31,  Whitened  station  and  network- summed  spectra 
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Figure  32.  Unwhitened  station  and  network- summed  spectra 
KNICKERBOCKER. 


59 


NETWORK 


'■k. 


SUM  0 


T - 1 - i - I  I  T 


FREQUENCY  (HERTZ) 


Figure  33.  Whitened  station  and  network- sununed  spectra 
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Figure  35.  Whitened  station 


and  network-summed  spectra 
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Figure  36.  Depth  phase 
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APPENDIX 

EARTHQUAt'c  SPECTRA 
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An  analysis  of  the  network- summed  spectra  for  17 
earthquakes  recorded  at  telcseismic  distances  (Tabl  I) 
indicates  that  the  spectra  for  six  events  exhibited 
spectral  nulls  which  may  be  related  to  interference  of 
the  type  pP-P,  Examples  of  these  spectra  are  given  in 
Figure  1.  Instrument  response  was  removed  from  the 
individual  station  spectra  prior  to  network  averaging. 

If  we  interpret  the  null  at  1.7  liz  in  the  spectrum 
of  the  16  March  1964  Kurile  Inland  event  as  being  the 
first  depth  phase  null,  the  data  imply  a  pP  delay  time 
of  0.6  seconds.  Assuming  an  over-burden  velocity  of 
3  km/sec,  this  delay  time  then  yields  a  source  depth 
estimate  of  0.9  km.  A  similar  argument  applied  to  the 
broad  band  null  at  2.5  Hz  in  the  spectrum  for  the 
25  April  1964  Kuriles  event  suggests  a  source  depth  of 
0.6  km.  Neither  can  be  considered  correct  in  view  of  the 
depth  estimates  of  140  and  40  km,  respectively,  reported 
for  these  events. 

Spectra  shown  for  the  Andreanof  Island  and  Chilean 
events  in  Figure  1  are  typical  of  the  network-summed 
earthquake  spectra  found  in  this  study.  For  completeness 
however,  network  summed  spectra  for  the  14  remaining 
events  studied  are  shown  in  Figures  2  through  4. 
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Figure  1.  Earthquake  spectra  network  sums. 
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Figure  2.  Earthquake  spectra  network  sums. 
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Figure  3.  Earthquake  spectra  network  sunx; 
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Figure  4.  Earthquake  spectra  network  sums. 
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